Applicability of Polytropes 


The four equations of stellar structure divide naturally into two 
groups: one describing the mechanical structure of the star us- 
ing (2.1.3) and (2.2.9), and the other giving the thermal struc- 
ture via (2.3.3) and (2.4.4). However, the only contact between 
the mechanical variables and thermal equations is through the 
temperature dependence of the equation of state. Under certain 
circumstances, however, the pressure can become independent of 
temperature, and only depend on density, t.e., 








P = Kø = Kpt!” (15.1.1) 


where n is the polytropic index. When the equation of state 
can be written in this form, the calculations of stellar structure 
simplify enormously. 


Decoupling occurs when: 


e Stars become supported by electron degeneracy. When this 
occurs, the pressure is related to the density by 
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(K = 1.2435 x 10°? i cgs) depending on whether the Fermi 
momentum is relativistic or not. For future reference, note here 
that the constant K depends only on u and atomic physics. 





e Energy transport comes only from convection, and radiation 
pressure is negligible. In this case, P ~ Pyas, and 
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When combined with the definition of Vaa, 
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this yields 
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or 
P = KpV/(-Vaa) (15.1.2) 


where K is not a function of atomic physics, but instead depends 
on the star’s boundary conditions. Note that for a completely 
ionized gas, Vaq = 2/5, which recovers P œ pre 


e Regions in which the ratio of the gas pressure to the radiation 
pressure is constant. If we let G = Pyas/P, then for an ideal gas 
plus radiation pressure equation of state 
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If we then use (15.1.4) to substitute for temperature, then equa- 
tion (15.1.3) becomes 
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where the value of K depends on the importance of radiation 
pressure. 








This last condition comes about naturally if the ratio KL/M is 
constant. To see this, consider that if the star is convective, P œ 
p!/ 0 -Vaa) as above, but if the star is radiative, then V = Vrad, 


and 
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If we then substitute using 


Vie (3.1.6) 
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Thus, if k£/M is constant, and the surface pressure is ~ 0, then 








(15.1.6) 
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This is the condition for the Eddington “standard” model. 


(Actually, setting «£/M constant is not a terrible assumption. 
Kramer law opacities decrease with temperature as T’~*-°, while 
the proton-proton chain gives en = L/M œx T°? for low mass 
stars. Thus, k£/M is approximately constant, and (15.1.7) is 
an adequate approximation. ) 





Polytropic Calculations 


To calculate the structure of a polytropic star, begin by assuming 
hydrostatic equilibrium, and multiplying both sides by r?/p, i.e., 
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If we then take the derivative (with respect to r) of both sides, 
and divide by r?, then we get 
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(Note that this is nothing more than the Poisson equation — if 
you substitute for pressure using 
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then you recover V?® = 4rGp in spherical coordinates.) If you 
then substitute in the polytropic relation (15.1.1), then 
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and Poisson’s equation becomes 
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To simplify this expression, let’s put it in dimensionless form. 
First, define a variable 6, such that 
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Poisson’s equation then becomes 
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Then define a dimensionless radius, € = r/r, where 
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where we have used (15.1.1) to substitute in the central pressure. 
In this form, the Poisson equation becomes 
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This is called the Lane-Emden equation; its solution gives the 
run of dimensionless density 0 as a function of the dimensionless 
radius, €. Since it is a second order differential equation, you 
need two boundary conditions. The first is at the center: from 
spherical symmetry, the pressure gradient at the center (0 = 1) 
must be zero. The second condition comes from the surface, £1, 
where the density should go to zero. So our boundary conditions 
are 


d0 


gE =0)=0 (the center) 


and 
6(€ =&,)=0 (the surface) 


Unfortunately, the Lane-Emden equation does not have an an- 
alytic solution for arbitrary values of n. In fact, there are only 
four analytic solutions. The first is for n = 0, which implies 
p(r) = pe, or a constant density sphere. For these models, 
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with & = V6 to satify the boundary condition of 0(€,) = 0. (This 
can be trivially checked via substition.) The second solution is 
for the n = 1 case, where 
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Note here that there are an infinite number of values of é; for 
which 0(£&1) = 0. However, in practice, £; = 7, since all other 
values would give an unrealistic zero density somewhere in the 
middle of the star. A third solution exists for n = 5, 
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For this model, 0 only goes to zero when £; —> ov, thus stars 
with n = 5 have infinite radius. All models with n > 5 have both 
infinite radius and infinite mass, but the n = œ solution is note- 
worthy. From (15.1.1), n = oo corresponds to P = Kp, which 
is the isothermal case. If we go back to the Poisson equation 
(15.2.2) and let n = oo, then the equation becomes 
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If we let Y = — ln p and £ = (4rG/K)!/?, then this translates to 
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This is the equation for an isothermal sphere. It has many ap- 


plications in astrophysics (especially extragalactic astrophysics). 


The remaining solutions of the Lane-Emden equation, including 
that for n = 3/2 (y = 5/3) and n = 3 (y = 4/3), must be com- 
puted numerically. For the inner part (€ < 1) of the polytrope, 
this can be done by writing 0 out as a power law, t.e., 
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and substituting into the Lane-Emden equation. When this is 
done, the left side of the equation becomes 


f(E) = 2a€7! + 6b + 12cé + 20d? +... (i5217) 


Meanwhile, the right side of the equation can be evaluated by 
expanding 0” as a Taylor series about € = 0. In other words, 
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F(E) = 0” = f(0) + f'(0)E + 
where, through (15.2.10), 
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The coefficients of (15.2.10) can then be found by equating the 
terms of (15.2.11) and (15.2.12). In other words, since (15.2.12) 
does not contain a €~! term, a = 0, and 
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Note that all the odd power terms are zero, due to the fact that 
a = 0. This leaves 
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The outer regions of the polytrope can now be computed numer- 
ically by translating the second-order Lane-Emden equation into 
two first-order equations with known starting point boundary 
conditions. In other words, if we write Lane-Emden as 
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dg dg £ 
This can be solved in a straightforward manner by starting from 
the polynomial expansion and using Runge-Kutta integration. 


then 


Some key values resulting from this integration are tabulated 


below. 
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